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We study the double scaling asymptotic limit for large matrix dimension of 
the partition function of the unitary ensemble with weight 



w{x) := exp --^ H - 

We derive the asymptotics of the partition function when z and t are of 0{N^2). The 
■ results are obtained using the Deift-Zhou steepest descent method and are expressed 

, in terms of a solution of a fourth order nonlinear differential equation. We also 

compute the asymptotic limit of such a solution when zA^2 0. The behavior of 
this solution, together with fact that the partition function is an odd function in the 
variable t, allows us to reduce such a fourth order differential equation into a second 
' order nonlinear ODE. 
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§ ; 1 Introduction 



In a previous article [23j we studied the asymptotics of the partition function 



^^^"'')^=MLn-p(-^+^-y) n w-^^'d^^ (1.1) 

i=l ^ 3 ^ ^ l<j<k<N 

in the range ciN~^ < z < C2N^, where ci and C2 are two constants independent of N. In 
this paper we extend the results in [23] to cover the case when the parameters t and z are 
of the following order 

t = -^(l + 0(iV-^)), ^= W^(l + 0(iV-|)), asiV^oo, (1.2) 



*The authors acknowledge financial support by the EPSRC grant EP/G019843/1. 
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where u\ and U2 remain finite. 
Partition functions of the form 



I r ^ 

iVT n^*^^i) n W-x,j\^d^x, (1.3) 



iN 



where w{x) > is the weight and J C M or J = S^, are ubiquitous in Random Matrix 
Theory because they contain all the information on the correlations of the eigenvalues. 
The asymptotic properties of the integral (11. 3p are determined by the analytic properties 
of w{x)] the weight of the partition function (11.11) has an essential singularity at the origin, 
which makes its computation particularly difficult. 

Besides being an interesting mathematical problem in its own right, the knowledge of 
E]sf{z,t) is important for several reasons. It determines the value distribution of 



where 

N 

An{x) ■.= Yl{x-Xj) (1.5) 

i=i 

and Xi, . . . ,xn are real random variables governed by a given joint probability density 
function (j.p.d.f.). Berry and Shukla [1] discovered that the value distribution of Qa? is a 
sensitive measure of the rigidity of the zeros of A^{x). Furthermore, they showed that all 
the moments of the distribution of can be computed from (II. ip . 

The function (II. 4p (or more precisely the denominator of Qn) was introduced by 
Tuck [25] in a numerical study of the Riemann zeta function. Tuck replaced Aat by the 
Hardy function 

Z{t) := t^^exp (^) n-h^r (i + i^t) C + ^t) , (1.6) 

where ({z) is the Riemann zeta function. He realized that the Riemann hypothesis implies 
Qn > 0. (Note that the opposite is not true.) 

The quantity EN{z,t) can also be used for other purposes. For example, if we set 

s = z^l2 = -t, s> 0, (1.7) 

then N\En{\/2s, —s) becomes the moment generating function of the singular linear statis- 
tics 



fit Wj 



The function 
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can be thought of as a singular perturbation of the GUE weight e~~, in which an essential 
singularity emerges. The asymptotics of the partition function with the weight f ll.9p can 
be computed using the method of orthogonal polynomials. Then, one needs to study the 
asymptotics of the polynomials orthogonal with respect to the weight (11. 9p . This can be 
achieved by using the steepest descent method developed by Deift et al in [15], [16] (see 
also [3). 

When the parameter z, and hence the singular perturbation, is large, the asymptotics 
of the orthogonal polynomials can be found with the help of the equilibrium measure, 
which is supported on two disjoint intervals. This was done in [23] (see also [9] .110]). 
The asymptotics of the orthogonal polynomials can then be used to find the logarithmic 
derivatives of the partition function with respect to the parameters z and t. However, as 
the parameters z and t become small, the gap between the two intervals in the support 
of the equilibrium measure becomes small as the singularity disappears. Therefore, the 
results in [23] no longer provide an accurate description to the asymptotic behavior of the 
partition function (11.11) . The purpose of the current paper is to obtain the asymptotics of 
Ei\i{z,t) when the parameters z and t are small. This requires the use of Painleve type 
nonlinear differential equations. 

In order to obtain the asymptotics of the orthogonal polynomials with weight (11.91) . 
in ^we express the orthogonal polynomials as a solution of a Riemann-Hilbert problem. 
In §1] the equilibrium measure for the GUE is used to modify such a Riemann-Hilbert 
problem into explicitly solvable 'model Riemann-Hilbert problems' in different regions of 
the complex plane, which allow us to study the asymptotic behavior of their solutions. 
Since an essential singularity emerges at the origin, the model Riemann-Hilbert problem 
in a small neighborhood of the origin is solved using a local parametrix constructed from a 
special solution of a nonlinear differential equation (^. This technique is widely used in 
the studies of double scahng limits in Random Matrix Theory (see e.g. [8], [11], [T2], [13]). 

Usually, the local parametrix is expressed in terms of a special solution of a nonlinear 
differential equation, which in turn gives the compatibility conditions for the following 
system of linear ODEs: 

d^^Cu) = A{C,uMC,u), 5„$(C,m) = BiCuMCu) = 0, (1.10) 

where m is a parameter measuring the strength of the perturbation from a critical point. 
In our case the parameter u is replaced by Ui and U2 in (II. 2p . The matrices A{() and B{() 
are rational functions of (. 

The compatibility conditions for (ll.lOp are among isomonodromic deformations which 
include the Painleve equations. In most of the previous studies of double scaling limits 
appearing in Random Matrix Theory, one parameter was sufficient to represent the whole 
family of perturbations from the critical point. Hence the conditions in (ll.lOp gave rise to 
a nonlinear ODE describing the double scaling limit. In most cases such double scaling 
limits are described by the Painleve equations. 

In our study, however, the local parametrix turns out to be more complicated. Indeed 
compatibility conditions analogous to (ll.lOp give us a nonlinear PDE instead. By using 
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the Hamiltonian theory of isomonodromic deformations, we are able to express such a 
PDE in terms of nonhnear ODEs in the variables ui and U2- This simplifies the problem 
considerably from the computational point of view. The ODEs that we obtain are of fourth 
order, and have not been appeared in the studies of Random Matrix Theory before. The 
Hamiltonian approach used in this paper can also be adopted in other more general cases 
of double scaling limits to express local parametrices in terms of solutions of nonlinear 
ODEs (see Theorem lO) . 

After deriving the nonlinear ODEs describing the local parametrix, in ^ we study its 
behavior as M2 — ?■ to obtain the initial conditions. Such initial conditions, together with 
the fact that the partition function E]^{z,t) in f ll.ip is odd in t, allow us to simplify the 
nonlinear ODEs further and describe the series expansion of the partition function in t by 
a solution of a second order nonlinear ODE (^ ITO|) . Since this ODE is a reduction from 
an isomonodromic deformation that satisfies the Painleve property, it should satisfy the 
Painleve property too. Hence, it should be equivalent to one of the Painleve equations. 
However, we have been unable to transform it into a Painleve equation. 

2 Statement of results 

The average f ll.ip is an entire function of t, therefore its Taylor series has an infinite radius 
of convergence and we can write 



In this paper we will assume the parameters z and t are of the order indicated in fll.2p 
and compute the leading asymptotics of the coefficients E^rniz) in f l2.ip . The case when 
z is of larger order was considered in [23]. When z and t are in the range studied in 
this article, the ensemble average f ll.ip is given in terms of a solution of a fourth order 
nonlinear differential equation. Such an equation can be expressed in terms of a time- 
dependent Hamiltonian system and is related to a particular isomonodromic deformation. 
The Hamiltonian equations are given by the following. 



00 




(2.1) 



m=0 
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Theorem 2.1. Let ui and U2 he given by and Hi, H2, hi and /i2 by 

Hi = 2i— + \u1Q2 - I-U2Q1Q2 - -^ujQi + \uiQl - 

U2 2 2 Au2 2 4 

+ -P1P2Q2 + -PIQ1, 

U2 U2 

H2 = + + -U2QIQ2 + -U2QI - —PIQ2 

U2 U2 O 8 2U2 



2U2 U2 8U2 4m2 



(2.2) 



8 -Ug 2^2 8^2 

, _Pi , _ P2Q2 PiQi ui ^ 

m — — , ii'2 — 773-^1' 

U2 U2 2U2 2U2 



where the variables Pi, Qi, P2 and Q2 satisfy the following time- dependent Hamiltonian 
equations, 

Pip, rin.-Lh. Fin, pin.-Lh. 

J, A; = 1,2. (2.3) 



dPk _ dHj + hj dQk _ dHj + hj 



duj dQk ' duj dPk 

Then the logarithmic derivatives of the average in U.l\) are given by Hi and H2 as N ^ 00. 

«i F (2-4) 

oz 

The solution of the Hamihonian equations (12. 3 p can be characterized by a Riemann- 
Hilbert problem which in turn gives the behavior of the canonical variables Pi, Qi, P2 and 
Q2 as U2 — )■ 0. This is computed in §21 and the asymptotic behavior is given in Theorem 

It is often useful to write En{z, t) as a power series expansion in t. For this purpose, it 
is sufficient to consider the equation (12. 3p for j = 2 at mi = 0. In this case, the asymptotic 
behavior of the canonical variables Pi, Qi, P2 and Q2 is much simpler 



^ + 3 4 /2 ^ 32 ,,,,, 



4 (21 2I 3^2 2i\ i /16 64 



'5l=M A/— -2+3+ |— + — + -— M2 +0(m|;. 



P2 = Oiv" 



2)-. 

25 „ 4 /3-2t 7-2i 2i 2i \ i / 32 8 \ ^3 



^2 = ^ + 2 - - + I + — ^ - - _ U| + ^ N2 + 0(^2 



(2.5) 
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From (11. ip . it is easy to see that En{z, t) is an odd function in t. This symmetry condition, 
together with the relation between Hi and the logarithmic derivative of EN{z,t), allow us 
to obtain algebraic relations between the canonical variables Pi, P2, Qi and Q2. 

P2 = 0, Q2 = ^ . (2.6) 

As a result, when ui = 0, the four ODEs (12. 3p can be reduced to two first order ODEs. 
Theorem 2.2. Let Hq and ho be the following 

_ {ulQ^ - AP^f iP , _ PQ 

IIQ o 7 t^O • 

32^2 U2 2u2 

Then as N 00, the average U.l\) is given by 

E,(z,t) = (2.)f n ,! exp ( (- H,iu, - g (l + O(N-i) 

where the variables P, Q satisfy the following time- dependent Hamiltonian equation in U2, 

dP _ d{Ho + hp) dQ _ djHo + hp) 
du2 dQ ' du2 dP 

which is the following 

_ Q{ulQ^ - 4P^) P 

8^2 2U2 

^ PjulQ'-AP') Q I 

2U2 U2 

The canonical variables P, Q behaves as follows as M2 — ?■ 0. 





2M2 ^ 


-1 




TT 








4 


/ ' 


-2 + 




V M2vr 








The derivatives '^^ ]^^^ in the asymptotic formula of Ej\i(z,t) can be expressed in terms 

of the canonical variables P, Q by using (12.31) and the substitution P = Pi, Q = Qi, 
together with the algebraic relations (12. 6p . For example, we have 

dui '"1=° ~ 4 ' 



This then gives the leading order asymptotics for the average (II. ip . 
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3 Preliminaries 

Let us introduce the variables Ui^n and U2^n by 

Ui,N = VNt, U2,N = Nz'^ (3.r 

The weight ( 11.91) then becomes 



WNiy) := exp| -A^ 

We can also define 



i/_„jge.p(-.(5^.|).^)j^J.-./.^, ,3.3, 

which is proportional to EN{z,t), namely 

En{z, t) = AT— Giv(Mi,7v, U2,n)- (3.4) 
Let us denote the poljTiomials orthogonal with respect to wiy{y) by TTj{y) and write 

N-l 

Gn{ui^n,U2,n) = det (/Uj+fe)^"^^ = J]^ hj, (3.5) 

j=0 

where 

/oo 
WN{y)y^dx, j = 0,l,... (3.6) 
-00 

and 

/oo 
WN{y)TTj{y)'r^k{y)dy = hjSjk. (3.7) 
-00 

Let us define the matrix valued function 

\KN-lTTN-l[y) ^^J_oo 

where kn-i = —2TiihM-i- It solves the following R-H problem: 
1. y{y) is analytic in C/M, 

3. i'fe) = (/ + 0(i/-')) (^p" !/^oo 

where 1 + (|/) and l^-(y) denotes the hmiting values of Yiy) as it approaches the left and 
right hand side of the real axis. It turns out that the partition function Giv(f i,f2) can be 
expressed in terms of Y{ii). 



8 



F. Mezzadh and M. Y. Mo 



Lemma 3.1 (Bertola, Eynard and Hanard [6j). (See also J2B^) The following differential 
identities hold: 



d\ogG 



N 



L-/ -Ti{Y-\y)Y'{y)a,)dy, (3.10a) 



dui^N AniN 

^^ = At-J -,Tr {Y-\y)Y\y)a,)dy, (3.10b) 
du2,N SmN^ Jy^Q y^ 

where the contour of integration is a small loop around y = oriented counter-clockwise. 

Note that, although the function Y{y) has a jump discontinuity on the real axis, it has 
a uniform asymptotic expansion near the origin and the residues in the right hand sides of 
(13.101) are computed using this uniform asymptotic expansion. 

4 Riemann-Hilbert analysis 

In this Section we will apply the Deift-Zhou steepest descent analysis to the Riemann- 
Hilbert problem (13. 9p . 

4.1 First transformation of the Riemann-Hilbert problem 

Let us now define the function g{y) to be 

9{y)= / ^og{y - s)p{s)ds, (4.1) 

2 

where p{s) is the equilibrium measure for the potential V{y) = y'- 



Pis) = —VA^. (4.2) 



That is, we have 



2/ \og\y - s\p{s)ds - ^ = I, ye [-2,2], 
2/ \og\y-s\p{s)ds-\<l, yGM\ [-2,2]. 
for some constant /. Now let T{%j), g{y) and F{y) be 



(4.3) 



T{y) = e^r(y)e-^^(^)'^«e^, g{y) = ^ - g{y), 

F(y) = -^( - ^ 
2 \2N^y^ Ny 



(4.4) 
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Figure 1: The opening of the lens around the interval Sj, j 
are the regions between 'E±j and the real axis. 



1,2. The lens regions L±^ 



Then T{y) is a solution to the following Riemann-Hilbert problem: 

1. T{y) is analytic in C \ M, 

2. T+{y)=T^{y)My), yeR, 

3. T{y) = I + 0{y-'), y ^ oo, 

where Jrin) is the following matrix on M. 



My) 



--N(g+{y)-g.iy)) ^-N {g+{y)+g.{y))+2F(y) 







,N{g+{y)-g-{y))) 



y e 



4.2 Opening of lens 

We now perform a standard technique in the steepest decent method (see [7|, [U], [IS])- 
Let yo be a point yo = 0{N~^) and define the lens contours around [-2,2] as shown in 
Figure [TJ Let the matrix S{y) be 



r ny), 



s{y) 



y outside the lens. 



^iy) [ ^2Ng(y)-2F(y) i ) ) ?/ ^ 



_Q2Ng{y)-2F{y) ^ j > 

1 0' 



(4.5) 



^2Ng(y)-2F{y) i I 5 



yeL. 



for j = 1,2. Then, the function S{y) satisfies the following Riemann-Hilbert problem: 



1. Siy) is analytic in C \ R; 

2. S+{y) = S^{y)Js{y)^ y G M; 

3. S{y) = I + 0{y-^), y ^ oo. 



(4.6) 
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where Js{y) is now defined by the following 

1 0\ 

^2Ng(y)-2F(y) ^ I ' 



Js{y) — [ ^2Ng(v)-2F(v) 1)5 V ^ ^±,ji J " 1)2, 



"■2.N , "l.iV - 



e 2^^a^ '^J' \ 

Js{y) = I "2,iv , ye [-2, -yo] U [yo,2], (4.7) 

g 2]v2y2 iVi/ / 

piV(9+(s/)-9-(j/)) p-2iV3(y)-2F(y) \ 



By the inequalities (14. 3p . we see that outside of some discs D±2 and Dq of radius 5 centered 
at the points ±2 and zero, the jump matrix Js{y) can be approximated by 



"2, AT , "l.JV ■ 



Js(y) = .J ^ ^""^ "• =(/+ o (iv-)) ( « ;y (4.8) 

This suggests the following approximation to S{y) outside of Dq and D±2: 

1. S°^{y) is analytic in C \ [-2, 2]; 

2. 5-(y) = 5-(y)(_°^ J), yG[-2,2]; (4.9) 

3. S^iy) = I + 0{y-'), y ^ oo. 

The outer parametrix that solves (14. 9 p can now be constructed as in [I5j and [H] : 



1 

7+7 " 7-7 ' \ /y_2\4 



5°°(y) = J^-. , 7 = ^ , (4.10) 



7—7 7+7 

2i 



,2/ + 2, 

where the branch cut is chosen to be on [—2, 2] and 7 ~ 1 as y — t- cxd. 



5 Local parametrix near the origin 

Near the origin, the approximation of S{y) by S°°{y) fails and we must find an exact 
solution to the Riemann-Hilbert problem for S{y) and match it with S°°{y) as y moves 
away from the origin. Let S be the union of the lens contours, then inside the neighborhood 
Dq, we want to solve the following Riemann-Hilbert problem: 

1. S^°\y) is analytic in Dq \ {Dq n (M U S)), 

2. sf{y) = S^'\y)Js{y), y G Do n (M U S) , (5.1) 

3. S^'\y)={l + 0{N-'))S^{y), z e dD,. 
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where the jump matrices are given by 

Jsiy) = { mg(y)+^-^ yeE^.nDo, j = i,2, 

Js{y) = { "2 iv "i.iv , y e {-6,-yo)U {yo,6), (5.2) 

_g 2]V2y2 ATy / 

"2, AT I "l.AT \ 



5.1 Conformal maps inside Dq 

Let ( be the analytic continuation of the following map in Dq: 

C = -iN {g+{y) - ^+(0)) = N{y + 0{y^)) . (5.3) 

Note that ^+(0) G ?R. 

For small enough 5, C is a conformal map inside Dq that maps the boundary of Dq into 
the point C, = oo a.s N oo. Let us define U2{y) and Ui{y) in Dq to be 

^2(1/) = M2,7V^^^ = M2,7V(1 + 0(?/^)), tti(?/) = = Mi_^(l + 0(?/^)). (5.4) 

Then, for small enough 5, both these maps are conformal inside Dq and we have 

^^^^ ^ 2iVV ~ ^ " —±^Ng+iO), ±Im(C)>0. (5.5) 

Let r be the union of the contours in the right hand side of Figure [2J Consider the following 
Riemann-Hilbert problem in the ( plane: 

1. -P(C) is analytic in C \ F, 

2. P+(C) = P-(C)^p(C), CeF, 

3. PiC)={l + OiC'))(^\ J), Im(C)>0, C^oo. ^^-^^ 
P(C) = / + 0(r^), Im(C)<0, C^oo. 
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Figure 2: The map C, that maps the neighborhood Dq into the complex plane. 



where the jump matrices Jp{C) are given by the following 



MC) 
MO 
MO 
MO 



0^ 



1 0' 



e ^ « 

"2(i/) "1(a) 

.-e c 

"2(a) I "i(a) ' 



C e U Te, 

Cerour4, 



(5.7) 



Then the local parametrix in (15. ip can be obtained from the solution of (15. 6p by 

^(°)(y) = 5°^(y)e-^^+(°)'^^ (J "q^) P(C)e^^+(°)'^^ lm(y)>0, 

as Ng+{0) G ^M. 

5.2 Existence of the local parametrix 



(5. 



We will now show that the Riemann-Hilbert problem (15. 6p has a unique solution. First let 
us define 

[ P(0, Im(C) < 0. 

Then -P(C) is a solution of the Riemann-Hilbert problem 



1. -P(C) is analytic in C \ F, 

2. p^(c) = P-iOMO, CeF, 

3. P{0 = i + O{c'), C^oo, 



(5.9) 
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(5.10) 



where the jump matrices JpiC) are given by 

( "2(h) I \ 
iiiM-iiiM 1 ' C ^ To u 
« y 

Since the Riemann-Hilbert problem and the jump matrices have no singularities at the 
self-intersection points ±yo of the contours, and since the jump matrices satisfy the cyclic 
conditions 

at the points ±?/o, where Jp^ is the jump matrix on F^, the theory in can be applied 
to (15.91) . In particular, (15. 9 p has a unique solution if there is no non-trivial function -Po(C) 
that satisfies 1. and 2. in (15. 9p and behaves as -Po(C) = 0{C,~^) as C ^ c>o. Suppose that 
such a function exists and let the region between Fi (F3) and Fq (F4) be f2o,+ (^4,+) and 
the region between Fq (F4) and Fg (F5) be f2o,- (^4,-)- Then, define the function X{Q to 
be 

^(0 = ^(0, Cec\^], 

J y ' j, C e (^0,+ u ^4,+, 

. I 1 o\ 

X(C) = Po(C) .2ic+^-^ J' C e ^^o,- u ^4,-, 

\ — e ^ c 1 y 

where = fio,+ U fio,- U ^4,+ U ^4,-. Then as in [16j, X{Q has the following properties: 
1. X{Q is analytic in C \ M and continuous down to M, 



3. X(C) = 0(r^), C^oo 



The function X{C)X{C,) is analytic in C+ and of order 0{C, ^) as C ^ 00. Therefore, by 
Cauchy's theorem, we have 
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Figure 3: The regions in the C, plane. 



By adding this matrix to its Hermitian conjugate, we see that 



'"2(«) I 



This imphes that the first column of X^{Q is identically zero. From the jump conditions 
in (15. lip , we see that the second column of X^{C) is also identically zero. Now let the first 
column of X be Xi and the second column of X be X2 and define the vector f{C) by 



/(C) 



^2(0, 



Im(C) < 0; 
Im(C) > 0. 



Then from (15.111) . we see that f{Q is analytic in the whole complex plane and behaves as 
as C — 00. Hence by Liouville's theorem, we have /(C) = 0. This shows that there 
is no non-trivial function Pq that satisfies 1. and 2. of (15.91) such that Pq = 0(C~^) as 
C — 00. Therefore, the Riemann-Hilbert problem (15. 6p is uniquely solvable. 



5.3 Painleve type differential equations 



We will now transform the Riemann-Hilbert problem (15. 6p into a Riemann-Hilbert problem 
with constant jumps but with essential singularities at C = and C = 00. The theory of 
isomonodromy deformations developed by Jimbo, Miwa, and Ueno can then be applied to 
derive differential equations that give the solution to (15. 6p . 

Firstly, let us deform the jump contours of (15.61) as in Figure [31 Let the solutions of 
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fl5.6p in the region I be Pi{C), etc., and define the function $(C) in different regions to be 

-e + c \ 

^ , (eluiv 

c / 

/ _e ^c^+ c \ 

Ve^- c ; (5.12) 

$(C) = P(C)e-'^^^^ C e ///, 

= p(() , c e \^ u y/// u ra, 

Let F-t be the semi-circles in Figure [3] in the upper and lower half planes respectively, then 
one can check that the matrix $(C) satisfies the Riemann-Hilbert problem 

1. $(0 is analytic in C\(F+UF_ UF2), 

2. <l>+(C) = $-(C)-/<i.(C), CeF+UF_UF2, (5.13) 

3. $(C) = / + 0(r^), C^oo, 

where J<i>{C) s-re the matrices given by 

1 0\ 

Mx) = I _2ic+^-^ J . C e F_, 

1 e ^c^+^^ \ 
J$(a;) = (a) ^liv) 5 C ^ r+, (5.14) 

-e c y 

J^(a:)=|J ^ ^ ^ I, CeF2. 
The matrix \&(C) defined by 

^(0 = $(()er-^+— i"^ (5.15) 
satisfies a Riemann-Hilbert problem with constant jumps on F+ U F_ U F2: 



'^'I'(^) =(01)' ^ ^ ■'"2' 



(5.16) 
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and has essential singularities at the points C = and Q = oo, 

no = MOe^'^^—^''\ c^o, ^^^^^ 

where ^o(C) and ^oo(C) have the asymptotic expansions 

^o(C) = ^r + M°^C + ..., ^oo(C) = / + M°°V^' + ---- (5.18) 

5.4 Differential equations 

From (15 ■4p . we see that Uj{y) behaves as follows as — )■ oo: 

«.(?/)=«.,^(l + o(|^)) =^,(l + o(^)), j = l,2, 

where Uj are the variables in (11.21) . Let $c(C) be the solution of the Riemann-Hilbert 
problem (I5.14p with Uj{y) replaced by Uj. Since the jump contours of ^(C) are all bounded, 
the matrix $(C)$^^(C) will have jump discontinuities on the contours r2 and r± of order 
I + 0{N~'^) and behave as I + 0{(~^) as ( ^ oo. By the general theory of Riemann-Hilbert 
problems, (see e.g. [TB], §7 and also §7]) we have 



^(^^^^(C) = I + O ( I . 

c lU +^^^^2(1(1 + 1); 



Hence, the matrix $c(C) approximates $(C) up to 0(A^^^) as — t- oo. From now on, we 
shall replace the functions Uj{y) by the constants uj and denote the solution to (I5.14p and 
(I5.15P with Uj{y) = Uj by the same notation <l>(^) and ^E'(C)- 

Let us consider the Riemann-Hilbert problem with constants jumps (I5.16P and asymp- 
totic behavior (I5.17p . Since the matrix \E'(C) has constant jumps on the complex plane, the 
matrices defined by 

^'^^QvI;-i(0=A(C), ^^^-\C) = BjiC), J = 1,2 (5.19) 



dC ' du 



'J 



are rational functions with poles at C = and C = oo only. By using (15.171) . we see that 
they are given by 
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where App denotes the singular part of A at C = 0. 
We can write the matrix A as follows. 



^(0 



(° 'j)c-' + -3. (5.21) 



where the form of the coefficient of ( ^ can be seen by expanding the first equation of 
fl5.19l) near ( = oo, 

= ,»3 + !|*;"',»3ir' + o(c-'). 

From f l5.20p . we see that the coefficients of s-^id iii the determinant of A{() 

coincide with the coefficients in the determinant of | — 0-3. This implies 



al + hcs = ^uj, 20203 + 62C3 + C263 = -^UiU2, 



1 (5.22) 

biCs + C163 + 62C2 + al = -uf. 

The compatibility of the linear differential systems in fl5.19p implies 

d,A{O-dcB,{O + m),B,{O] = 0, J = 1,2, 

d,B,iC)-d2B,iC) + [B2iO,B,iC)] = 0, ' 

where dj is the derivative with respect to uj and the derivative with respect to (. 
Although these compatibility conditions can readily be expressed as PDEs in the variables 
Ml and U2, it will be much easier to study the variations if we can write down ODEs in 
these variables. 



5.5 Hamiltonian equations and ODE in the variables 

We will now express the compatibility conditions f l5.23p as a system of Hamiltonian ODEs. 
First let us define a Poisson bracket on the space of matrices of the form A{() = X]j=i ^jC~''+ 
ia^. Let the loop algebra q be the set of smooth maps / : S*^ — )■ 512(C). We will split q into 
subalgebras 0+ and 0_, where g+ (g_) is the set of maps that admits holomorphic extension 
to the inside (outside) of the unit circle and that the maps in g_ vanish at infinity. Then 
by using the pairing < X{(),Y{() >= Res^=o Tr (Xy), we can identify g_ with the dual 
of g+. This gives the Poisson brackets on g_ defined by (see, e.g. [T7] ) 



{f,g}{A) = {A, [df{A),dg{A)]), 



(5.24) 
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where / and g are functions on g_ =0+- By using the method of the moment map (see 
[I]; [3], [E] and [IH]), this Poisson structure can be restricted to the space of matrices of 
the form A[Q = Yl^j=o^j^~'' ~^ ^"^s- '^^^ Poisson brackets between the matrix entries aj, 
bj and Cj are given by 

{a ., bk} = -bj+k~i, {dj, Ck} = Cj+k-i, j + k <4:, 

(5. 25) 

{bj, Ck} = 2aj+k-i, j + k-l <A 
and all the other Poisson brackets vanish. 

Remark 1. Note that, although strictly speaking, the matrix A{() = X]j=o^iC~"' + ^crs 
does not belong to the Lie algebra 0-_ because of the ia^ term, one can think of the space 
of matrices A as being parameterized by the entries aj, bj and Cj and use (15.251) to define 
a Poisson bracket on this space. In fact, the original construction in [T] and [3] applies 
to much more general form of matrices and readers who are interested may consult these 
references. 

The symplectic leaves in this Poisson manifold are the co-adjoint orbits of the matrix 
A, specified by det{A) = constant. A set of canonical coordinates on these symplectic 
leaves can be found as follows (see [3] and [21]). Let i'iC) be an eigenvector of ^(C) with 
eigenvalue A(C). Then, the poles qj of ipiC) and the values X{qj) of the eigenvalue A at qj 
satisfy the relations 

{Pj, qk} = 5jk, {Pj,Pk} = {qj, qk} = 0. 
For example, let ip{C) be an eigenvector of A{() normalized by 

(1,0) -^(0 = 1. 

The poles qj of ipiC) are then the zeroes of the polynomial 

he + + &3 = 0, 

and the eigenvalues at these points are given by 

3 

Pj = - ^(^kqj'' -i- 

k=2 

Therefore, the symplectic form on these symplectic leaves are given by 

2 Wl _ _ \ 

u = — dakqj^ A dqj = ( -da^ A dqj'^ + da2 A dq~^ j 

j=l k=2 j=i ^ ^ 

b2 \ J b2\ , J bi 
= a —a^- — h 02 A rf — — + da^ A a 



\ 'h J \ bsj \ h. 

Hence, the following gives a set of canonical coordinates on these symplectic leaves 

^1 = -^37^ + 0.2, Qi = P2 = as, Q2 = -r"- (5.26) 
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A set of coordinates on the Poisson manifold is then given by -Pi, Qi, -P2, Q2, ^i, ^2- The 
deformations A = [B,A\ can be expressed as Hamiltonian flows in this Poisson structure 

0, m- 

Proposition 5.1. Let B = {(''A)_, where is the projection onto then the following 
equation 

A = [B, A] 

is Hamiltonian with respect to the Poisson bracket ^^572^. Its Hamiltonian is given by 

iJB = ^Resc=oTr(AV) 
In particular, the flows A = [Bj, A] have the following Hamiltonians. 

H, = -^Resc=o Tr {A'C^) , = --^Resc=o Tr (Ia\] - ^H,. (5.27) 
2u2 2u2 \2 J 2u2 

By using fl5.22p . we can eliminate the cj and express the Hamiltonians in terms of the 
canonical variables fl5.26p . They are given by (12. 2p . Due to the presence of the terms 
d(^Bj in first equations of (I5.23p . the deformations of the matrix A are not of the form 
A = [Bj,A\ and hence the Hamiltonian flows generated by (12.21) are not sufficient to 
describe the deformations of A. However, by using a technique developed in [22], we can 
write the first equations in (I5.23P as a time-dependent Hamiltonian system. First let us 
choose a set of coordinates i/i, . . . , z/4, -ui, ^2 on the Poisson manifold such that the partial 
derivatives di and 82 (derivatives that keep z/i, . . . , 1/4 fixed) satisfy the following equations: 

d^a2 = , d^b2 = , <9iC2 = , 

U2 U2 U2 





= dlb. 


= dlc^ = 0, 




= ^r&i 


= d\c^ = 0, 




_ cts 
M2' 


d^h = ^, 

U2 


dla2 


1 

~ 2^2 


( ^""^ \ 

a2 H as 

V U2 J 


0^202 


1 

~ 2^2 


( ^""^ \ 

C2 H C3 

V U2 



die, = ^ 



C3 

M2' (5.28) 



d!^b2 = — {b2 + -b, 
2U2 \ U2 

dl^ai = d^bi = d'^ci = 0. 

By using (I5.22p . one can check that this definition is compatible with djUi = 6ij. Note 
that the derivatives (I5.28P imply the following 

d'(A = d^Bi, d^A = d(;B2. (5.29) 

Therefore we can write the equations (15.231) as a time-dependent Hamiltonian system, 

djA = {Hj,A} + d^A. (5.30) 
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In order the find the evolution equations in the form djPk = ~qq^j djQk = -g^ for 
suitable Hamiltonians Hj, we must replace the partial differentiation dj by the partial 
differentiation dj""^ that keeps the canonical coordinates fixed and modify the Hamiltonians 
Hj accordingly. We have, by the chain rule, 

2 

QU ^ Qcan ^ ^ ^gup^g^^ ^ d^Q.dQ,) . (5.31) 
k=l 

From f l5.28p and f l5.26p . we obtain the derivatives d':^Pk and d'^Qk- 



dlP, = dlP2 = dlQ2 = 0, ^',Q^ = -, 

U2 



92^P2 = -, ^2 = - — , d!^P^ = ^P„ d^,Q, = -^Q,-^, 
From this and f l5.3ip . we have 



(5.32) 



, _Pi , _ P2Q2 PiQi p (5.33) 

f^l — — ) '^2 — TT^-^l- 

U2 U2 2^2 2-U2 

Therefore (15.30p becomes 

djA = {Hj + h,,A} + df''A. 
In particular, we have the following. 

Theorem 5.1. The deformations of the matrix A in Ii5.21\) 

d.AiO - d^B.iO + m), B,iC)] = 0, j = 1, 2, 

are equivalent to the time- dependent Hamiltonian equations 

d{H, + h,) d{H, + h,) 

d,Pu- d,Q, - , j,k -1,2 

where Hj and hj are given by h2. ^) and / (5. 5*31) . 

These give us ODEs in the variables Ui and U2. A special solution of these ODEs, 
specified by the jump conditions (15.161) . will give us the local parametrix via (I5.19p . (I5.15p . 
(I5.12P and (15. 8p . In particular, we will show that the leading terms in the logarithmic 
derivatives of the Hankel determinant Gn with respect to ui and U2 in (13. 3 p are given by 
the Hamiltonians —Hi and —7^2- 



6 Local parametrices near the points ±2 

Near the points ±2, the approximations of S{y) by S°°{y) also fail and we must again find 
exact solutions to the Riemann-Hilbert problem of S{y) and match them with S°°{y) as 
y moves away from these points. Such local parametrices can be constructed using Airy 
functions as in [S], [15], [IS] and we shall not repeat the construction here. 
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7 The final transformation of the Riemann-Hilbert 
problem 

We now show that the parametrices we constructed in §4.21 and §3 are indeed good approx- 
imations to the solution S{y) of the R-H problem fl4.6p . 
Let us define 

where S^^^y) are the local parametrices inside the neighborhoods D±2 and Dq. Then, the 
function R{y) has jump discontinuities on the contour shown in Figure |H In particular, 




Figure 4: The contour F^. 

R{y) satisfies the R-H problem 

1. R{y) is analytic in C \ Tji, 

2. R+{y) = R4y)My), yeTR (7.2) 

3. Riy) = I + Oiy''), y ^ oo. 

From the definition of R{y) it follows that the jump matrix July) has the following order 
of magnitude as — > oo: 

[I + Oie-^"), for some fixed r/>0 for ?/GFi?\((9D±2 U UFo), 

where Fq is the part of [—2, 2] that is outside of D±2 and Do- 
Then, using well established techniques (see,e.g' [I5] §7), we have 



= ' + (7.4) 

uniformly in C. Therefore, the solution S{y) of the R-H problem fl4.6p can be approximated 

by S°^{y) and S^P^y): 

f(/ + 0(iV-i))5(P)(y), yeD^^l^Do, 

Combining these expressions with equations (14.51) and (14. 4p we obtain an asymptotic for- 
mula for the solution of the original R-H problem (13. 9p . 
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8 Asymptotics of the Hankel determinant 

We will now use Lemma 13.11 to express the logarithmic derivatives of the Hankel determi- 
nant Gn in terms of the powers series $0 in f l5.17p . These derivatives will then be identified 
with the Hamiltonians Hi and H2 in f l2.2p . Firstly, from the definition fl4.4p and f l4.5p of 

T(?/) and Siy)^ we see that 

Y{y) = e^S{y)e('-^''^y^-^>\ (8.1) 

From dUD, we have S{y) = R{y)S^^^{y) near the origin, where R{y) a matrix of order 
indicated in (17. 4p . Then from (15. 8p we obtain 



y; = e^R{y)S'^iy)e-''<^+^">''-' ^ j p^(^)e^v^c;+(uj+g+(y)--j<x3^ j^^^^ > 

Y_ = e^i^(|/)^~(^/)e-^^+(°)"='P_(C)e^(^+(°)+^-(^)-^)"^ lm{y) < 0. 
Now from (15.120 and (15.150 . we obtain 



Therefore -P(C) are related to the expansion ^(C) by 

We can then use these to compute Tr (Y^^^^a^^ dy. 

Tr (^y-i^cTgj dy = 2Ng'^{y)dy ± 2idC + Tr {i>,\Qi>',{Qa,dc) 

+ Tr {k± (C) [S^ {y)Y\s^ {y))' K± (C) a,dy) 
+ Tr [k^\C) {S^{y)y' R'\y)R'{y)S^{y)K^{Ocr,dy 



From fia.lOap . fl3.10bp and the relat ion between y and C (15. 3p . we obtain 

= -^Resc.0^ Tr {^'{CWoiO^sdc) + 0{N-'), 
= ^Res,.o^ Tr (^0 ^(O^'olO^s^C) + 0{N-'). 

The leading order terms in these derivatives are the derivatives of the r-function for the 
equations (I5.23P |21| . Bertola, Harnad, Hurtubise and Pusztai have shown that they are 
given by the Hamiltonians (I5.27p . 
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Proposition 8.1. f2UV ) The Hamiltonians in l^5.21\j can he expressed as 

(8 3) 

Proof. First let A be the eigenvalue of the matrix A{C). Then the expansion of A near 
C = is given by 

A = ]^U2C^ - \u^C^ + ifi - + O(C'). (8.4) 
First note that, from fl5.19p . we have 

Let us now denote ^q^A^q by A and write it as 



2 VC^ Cv V^2i -n^ 

where 11 and Fij are bounded at C = 0. In particular, we have 

Tr {^,\C)m)^.) = p - ^ + 2n. (8.5) 

Let v± bet;± = A± ^|^|f — ^j+nj, then since A is an eigenvalue of A and A is 
conjugated to A, we have 

det (^A - = v+v. - F12F21 = 0. 

From the expansion of A in fl8.4p . we see that is of order 0(C~^), while -F12-F21 is of order 
0(1). Therefore, t;_ is of order C,^. Hence LI behaves as 

n = ifi-2if2C + 0(C') 

This, together with 08.51) . imply the proposition. □ 
This proposition and (18.21) implies (12. 4 p of Theorem 12.11 



9 Asymptotics of the canonical coordinates as 1^2 — ^ 



We will now compute the asymptotic expansion of the solution $(C) of the R-H problem 
(15.131) as M2 — )■ 0. We shall assume the scaling between the variables Ui and U2 is such that 
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ui = is finite. Let us first divide the ^-plane into 3 regions as in Figure [5l and define 
the function {Q to be 

<f(°)(C) = /, Cefil, 

$^°HC)=(j-.c t')' (9.1) 
Now let Aq be a small disc around the origin of fixed radius. Let 0(x) be 

1 r e-^+^i** 



Let ,^ = ^2 and define $(p)(C) to be 

*«(C)^p.,c W«-f-)^^-^), MC)>0. 



(9,3) 



where 0±j are the coefficients in the expansion of (j){x) as x — > cxd in the upper/lower half 
planes given by 

^ / s-^-^-^+^'^'ds. (9.4) 
'r± 



2Txi 



Here r+ (r_) consists of the two intervals (— oo, — e) U (e, oo), oriented from — oo to oo, 
together with the semi circle of radius e around the origin in the upper (lower) half plane. 
The exact value of e is not important as the integrand is analytic away from zero and 
hence by Cauchy's theorem, integrals with different values of e will yield the same result. 
In particular, by the residue theorem, we have 
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For the first few coefficients, this gives us 

0+,o-0-,o = l, 4>+,i- (t>-,i = Ui, 

Integrating by parts, we can express the coefficients in terms of 0±,o- 
Proposition 9.1. The coefficients (f)±j satisfy the following recursive relation. 

4>±,j = J (w~i0±j-i - 4>±,j-2) , j > 2. (9.7) 
The first two coefficients (f)±^i and (f)±^Q are given by 

(t)±,o = ±J + ^= / e^c^s, = ±^ + — ^ ( Ml / e^rfs - ) (9.8) 

^ V2vr Jo V 2vr \ Jo / 

Proof. The proof of (19. 7p follows from integration by parts immediately. 



s-i-ig-^+^i-cis = — I (ni - s) s-^e-^+"^^(is 

2m 2mj 



j ^ (Ml0±,j-1 - 0±J-2) 



for j > 2. For j = 1, the above gives 

(i)±,i = Ui(f)±^o + 7^-1 e~^+"'"c?s = Mi0±,o ^e^. (9.9) 

To compute 4'±fi, we expand the factor e"^* in a power series: 

k 



f -1 -slsr^ (uis) 



"^^± fc=0 

27r /r, 27r^ A;! 



r± 



k=l 



°° ~ 2fc+l 



2^yr/ ' ''^' + 7!^^;^ (2^ 1)2'= A;! 



fc=0 

ill 



2 



(9.10) 



= — / s e ^ ds -\ ^= / e 2 ds. 

27r Jo 

To determine the first integral in the above equation, note that if we change the integration 
variable from s to — s, then the path of integration will change from F-i- to — F^;. Therefore 
we have 

— / s-^e-^ds = / s-^e--ds. (9.11) 

2vr 271 
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Now, by a residue computation we see that 



— / s "'^ ds ^— I s "■^ ds = 1. 

27r /r, 27r 



Prom this and (19.111) . we see that ^ Jp^ s "2 ds = This, together with f lQ.lOp and 
(K9\f . imphes (EH]). □ 

From (19.61) . (19. 3p and the property of Cauchy transform, one can check that <^^p\C) 
satisfies the following jump condition 



"2 I "1 



$f(C) = $L^^(C) 

Hence, if we define -R#(C) to be the matrix 



then -R$(C) satisfies the R-H problem 

1. -R$(C) is analytic in C \ Tr^, 

2. /?$,+(C) = i?$,-(C)^R.(C), Cer^,, (9.13) 

3. /?$(C) = / + 0(r'), C^oo, 

where the contour Tpt,^ is indicated in Figure Eland the jump matrices are given by 

Jr.= \ ^ .2 M , C e Fh 

« 



1 0^ 



"2 / "2 I "1 



e ^ ^ (e ^ « - 1 

e-2< 1 - e ^+ O 2 - e 



(9.14) 



With the scaling ui = such that Ui is finite, the jump matrices Jr^ are of order 
/ + 0{^/u2) and hence, by the same reasoning used in §71 we have 



Asymptotics of a singularly perturbed GUE partition function 



27 




Figure 6: The contour Tr^. 



In fact, the matrix R<!,{C) can be solved as a series expansion in y/u^ by using the method 
in §7 of [IS]. In order to find the necessary initial conditions to solve the ODEs in Theorem 



15. 1[ we will need the expansion of -R#(C) up to the order m| . To do so, let us first simphfy 
the R-H problem of Rq, with the transformation 



RiC)=R^iC) 



i?(C) = i?$(C) 



q2(^ c — e^*^ ( e ^ 

_ '"2 "1 

e « 
1 0^ 
-2ic f e5c^"~ - 1 ) 1 



- 1 



, C e ^2 \ Ao, 



Then, the matrix R{C) satisfies the following R-H problem 

1. R{C) is analytic in C \ OAq, 

2. R^{C) = R4C)JRiC), CedAo 

3. R{0 = I + O{C'), C^oo. 



where the jump matrix is given by 



1 

2 "1"? 



E 



oo 

fc=l 



,fc Afc 



Eoo 



1 



E 



fc=l <p+,fc-^ 



fc 

2 



(9.15) 



where we have used (19.51) in obtaining the above equation. Note that is of order 
1 + 0(7^^). 

Let us denote the expansions of -R(C) and by 



fc=i 



fe=i 
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By comparing the coefficients of m| in the jump conditions of R, we obtain the following 
jump conditions for the coefficients in the expansion of -R(C)- 

= Jr.i + Ri-1 R+,2 = ( R-,iJr 1 + "^R 2 ) + 

This gives the linear jump conditions for the coefficients Rk- Together with the requirement 
that Rk = 0{C~^) as C — >^ oo, we can solve these R-H problems for the coefficients Rk- 
From (19.141) and (I9.6p , we see that ^ is given by 



Let £f be the following 



3 



£f{0 = C^-' - e^-^) = 0(1) in C. (9.17) 



.fc=0 



By using (19. 6p and (I9.16p . it is easy to check that the solution to R^^^ is 

Let us now denote the jump matrix of R)^ by 2 = R-,iJr,i + Jr,2- Then, we have 

J ^}_( 0-,i(0-,i-0+,ie-2<)+0_,2 02 ^ (e2< - 1) + 0_,2e2*^ 
^2 y^^^^ (0+,ie-2< - 0_,i) - 0+,2e-2< 0+,i0_,i (1 - e2<) - 0_,2. 

This gives 

_( -0+,i0-,i^r(C) + 0-,2) ^i+(c)\ . ^ . 
' " V('/>l,i - ^r(c) -0+,i0-,i^i'(c) J ' ^ ^ °' 

h <^-,2 - ^~10-,1 + 2z0+,i0„,iC 0-,2 + 22 (02 1 + 0__2) A > ^ . . 

(2 ^«~^0^_^_0^2 + 22(0+,2-0i,i)C -0-,2-2z0+,i0_,iC y*' 

(9.19) 



(9.18) 
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The entries of the jump matrix Jji^ = R-^iJf^2 + R-A'^R2 + ^^^s given by 



) 



where ao, a±^i and t±,o are given by 



'-zt.O = 20±,10±,2 =F 0±,3 - ^^10±,1- 



«±,1 = 0±,2 =F (/"I,! 



(9.20) 



Inside Aq the matrix R3 is 




ao^2-(C) -220+,ia_,i^r(C) i_,o^2+(C) + 2z0+,i02 i^+(C) 
6+,o^2"(C) + 2^0^,i0_,i^r(C) «o^2+(C) + 2z0_,ia+,i^+(C) 



) 



(9.21) 



and i?3 outside of Aq is given by the above plus 3. 

To compute the asymptotics of the coordinate Q2, we will need an extra term in R, 
which can be found in the same way as the other Rj: 



«±,4 = 0T,1 (0±,10±,2 =F 0±,3) - 0±,l(^T,O + "0) =F (y^,l4>±,2 

<-±,i = =F2(2a=p,i0±,i + ao)0=F,i, 

L±,2 = ±(pTA + (PT,2 (0=F,2 - «~10T,l) + 4>T,l i^pT.S ± (''T.O + "o)) , 

for C e Aq. 

We can now use this result and the expansion of to compute asymptotics of the 

matrix A{() in (15.201) as U2 — > 0. Firstly, note that from flS.lSp and (15.171) we have 



One can check that both definitions of in (19.31) give the same $0. From this and (15.201) . 
we can compute the asymptotics of the coefficients of A{Q. The computation is long and 
tedious but straightforward. They can be obtained via a computer algebra package such 
as MAPLE. The results are given by the following theorem. 



a±,2 




do = i?*(C)$^^HC)e^^' 
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Theorem 9.1. Let ui = and suppose ui exists and is finite as U2 — )■ 0. Then as 

U2 — )■ 0, the canonical variables Pi, Qi, P2 and Q2 in h2.2^) have the following behavior: 



p -^1 f , . , , A , / 20102 0i(0_,2-0_,i) 0_,l 
Pi = 77—^2 - ^ — ^ + 20-,o + 1 U2 + - —3 + -2 — 



V V-,0 / ^ 



+ 40_,i - 201 + 40_,i0_,o)^| + 2i(^ ^ ; + 



, 01 0-,1 (90'_,i + 9010-,! - 170_,2) - 0^0-.2 

,0 ' 903^0 



20„.i(30_,i + 0i) , 1001,1 - 50i0_,i + 0f , ^ ,2 , , 1 , ,^2 

' h D0_,i - 40i0_ ' ' "-^ 

1 / rh.rh . \ / 2010^;^ 01 ("^ — 6'^ 



+ qL' + o7 ^ + 601,1 - 40i0-,i + - + 40i,i0_,o U + 0(Mi 

901,0 30_,o 3 ' / 



5 



./ 20i02_i(190_,i + 170i) 401 (901 1 + 90102, 1 - 0i0_.2) 
, 60„,i (40_,i + 70i) , ^ 0i0,,i , 4 , 80_,o^ I 

H ■ T-T \- 8—7 + - H ^ )W2 + flMl) 



3 3 



ul 



P2 = (pU2 + i {4:(p^Ui -Ui+ 4(y90i) + (01^1 + (Sui^ - 1202) ^ _ 240^ii^<^2 _ X2ui2(^3^ 
Y.^2~ 8^~1^ /^A40Ui2-28\ ^^^3\ /I04ui=^ n.^2~^ 2 

+ t(^4<j)lui-^+ii j0, -320^^, + ( — ^ 960?Mij¥;2 

- 960iMi2y?3 - 32miV^)m| + 0{ul), 

^ o-<^-i --2 , . /^^^-.i , <^-,i-<^-.2\ ^./30i0l,i + 0l,i „010l,l , 2Mi0_,2 + 70_,i 
= 2«- — Ms +4 + 7 - 2t[ -3 + 6—2 — + — 

XI /0l 1 (3301 i + 790i0_,i + 3202) 602 ,(1102_^ + 140i0_,i + 2302)- 2020_,2 
+ ^2 - [-^ + ' 



-,0 



,1 



(220_,i + 1501 + 240i0l_i) 2401 1 + 520i0„,i . 5 lO0_,o 



+ 301, 30_,o ^ +^K. + -, + ^h + Oiul), 

where ip = f^^ e~^ds and the can be computed using Proposition \9.1\ In particu- 
lar, when Ul = 0, this gives ([H 



10 Reduction to a second order differential equation 

From fl3.3p . we see that the average G]^{ui^n^U2,n) is an odd function in ui^^ and hence 
the derivatives d'^'^^XogGN are zero when ui^n is zero. Then from (18.21) and (18. 3p . we see 
that 

52^ffi|„,=o = 0, A;g{0}UN. (10.1) 
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This leads to non-trivial algebraic relations in the canonical coordinates that cannot be 
easily seen from the Hamiltonian equations in Theorem 15.11 This allows us to express the 
power series expansion of Gat in ui^n in terms of the solution of fewer ODEs. 

From the first two equations in fllO.ll) . we obtain the two possible sets of relations 
between the variables 

P2 = 0, Q2 = - ^ (10.2) 

or 

iu\Q\ - %%ulPlQ\ + mPt - SnjPi ^ _ 4zPi 

The condition &l^Ji\\u^={) = does not yield any new constraints and the other equations in 
fllO.ip involve high powers of Pi and Qi which cannot be solved analytically. The relation 
(110.21) is compatible with the behavior of the coordinates Pi and Qi at U2 = fl2.5l) . while 
the relation (110. 3p is not. Therefore at ui = 0, the Hamiltonian equations in Theorem 15. II 
are reduced to the following 



Q{ulQ'-AP^) ^ P 

8U2 2U2 



(10.4) 



du2Q r, 3 o 

This can be represented as a set of time dependent Hamiltonian equations in the canonical 
variables P and Q with the Hamiltonian 

^_ julQ' - AP^ f iP PQ 

TL — T — — 

32^2 U2 2U2 

Since this system of Hamiltonian equations is a reduction from an integrable systems that 
satisfy the Painleve property, it should also satisfies the Painleve property and therefore 
be equivalent to one of the Painleve equations. However, we do not know how to transform 
it into a Painleve equations. 
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